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The influence of the rotation on the total masses and radii of the 
neutron stars are calculated by the Hartle's slow rotation formalism, while 
the equation of state is considered in a relativistic a — lv model. Com- 
paring with the observation, the calculating result shows that the double 
neutron star binaries are more like hyperon stars and the neutron stars of 
X-ray binaries are more like traditional neutron stars. As the changes of 
the mass and radius to a real neutron star caused by the rotation are very 
small comparing with the total mass and radius, one can see that Hartle's 
approximate method is rational to deal with the rotating neutron stars. If 
three property values: mass, radius and period are observed to the same 
neutron star, then the EOS of this neutron star could be decided entirely. 
PACS: 04.40.Dg 95.30.Sf 97.10.Kc 97.60.Jd 



I. INTRODUCTION 

Neutron stars are dense, neutron-packed remnants of massive stars after their super- 
nova explosions. Recently, in both experiment and theory, much concern is focused on 
the neutron stars [1-3]. One reason is to determine the equation of state (EOS) of super- 
dense matter, and then understand the early universe, its evolution to the present day 
and the various astrophysical phenomena; the other one is that neutron star is one of the 
more possible sources of detectable gravitational waves. To understand neutron stars, the 
first thing is to understand its structure, such as the compositions, total masses, radii, 
redshifts etc. Because of the strongly gravitational field and the high density, neutron 
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stars must be studied in tlie framework of general relativity. For a static neutron star, 
giving EOS, using TOV equations [4] , which are educed from Einstein field equations, 
one can get an exact solution. But for a rotating neutron star, the components of Einstein 
field equations become much difficult. Nowadays, several approximate solutions of this 
problem were developed [5,6], in which Hartle's slow rotation formalism [7,8] is the most 
popular one. 

The properties of neutron stars such as masses, rotational frequencies, radii, moments 
of inertia and redshifts are sensitive to the EOS of the matters [2,3]. As the interior core 
contains most of the mass of neutron star, so attention to the EOS is mostly focused on 
the core of neutron star, that is to the matters at density above several times of nuclear 
matter saturation density. In the core of neutron star, the compositions still keep blurry 
in some degree due to the high density and uncertainty of strong interaction, but as den- 
sity increasing in the neutron star, neutrons can drip out of nuclei and form a neutron 
gas; and due to the URCA process n ^ p + e~ + Pg and the neutrality of neutron star, 
there are at least protons and electrons in the chemical equilibrium; on the other hand, as 
the density increasing, hyperons will be dominant in the neutron stars [9]. In this letter, 
the Hartle's formalism will be used to deal with two kinds of rotating neutron stars: the 
traditional neutron stars, in which n,p, e, /i, are the main elements; and hyperon stars, in 
which n,p, e, A, S, S, A are the main elements. The EOSs of them will be considered 
in the relativistic a — uj model. 

In this letter, we adopt the metric signature - H — \ — h, G=c=l (the factors G and c 
will be reinserted correctly in numerical calculation). 



In relativity, the space-time geometry of a rotating star in equilibrium is described by 
a stationary and axisymmetric metric of the form 



where uj{r) is the angular velocity of the local inertial frame and is proportional to the 
star's rotational frequency fl, which is the uniform angular velocity of the star relative 
to an observer at infinity. Expanding the metric function through second order in fl, one 



II. HARTLE'S SLOW ROTATION FORMALISM 
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e^^ = 1 + - mo + m2P2 1 ^ " 1 ^ , 3 
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e^^ = r'sm'e[l + 2(t;2 - /i2)P2] , (4) 

e^''^r'[l + 2{v2-h2)P2\, (5) 

where e^"^ and Mo(r) denote the metric function and the mass of the nonrotaiting neu- 
tron star with the same central density, respectively; P2 is the Legendre ploynomial of 
order 2; the perturbation functions mo, m2, ho, /i2, ^2 are proportional to fl"^ and are to be 
calculated from Einstein field equations. 

From the {t, 0) component of Einstein field equations, one gets [7] 
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wh.eicej{r) — e~'^[l — 2Mo(r)/r]^, O — O, — lo, which denotes the angular velocity of the 
fluid relative to the local inertial frame. The boundary conditions are imposed as cu — cuc 
at the center, ^|^^=o , where u)c is chosen arbitrarily. Integrating eq.(6) outward from 
the center of the star, one would get the function uj{r). Outside the star, from eq.(6) one 
has uj{r) = Q — where J is the total angular mementum of the star, which takes the 
form [2] J — ^RQ^\r=Ro- So at the surface, one can determine the angular velocity 
corresponding to cUc as 

n = a;(i?o) + 2-^. (7) 

Prom the (t, t) and (r, r) components of Einstein fleld equations, one gets two coupled 
ordinary differential equations of mo and as [7,8] 

dmo ^ d{p + p) 1 .2 4.c^a) 2 1 ^ djf) _^ 

^ = -^^p + + 12^ ' - r ' 

dPo mo(l + 87rr^p) 47rr^(p + p) ^ 1 r'^j^ .did ^ 1 d^ r'f'uj^ ^ 



dr [r-2Mo(r)]2 r - 2Mo(r) 12 r - 2Mo{r)' dr ' 3drW -2Mo{r)' 

(9) 

where = — /io + |r^e~^"a;^-|-C , here C is a constant determined by the demand that ho 
be continuous across the star's surface. These equations are also integrated outward, with 
boundary conditions that both mo and vanish at the origin. With the same central 
density, the difference between the mass of the rotating star and the non-rotating star is 
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5M = mo(i?o) + ^, (10) 



The difference of tlie mean radius is 



5r^-pl{p + p)f^. (11) 



III. MODEL FOR THE EOSS - THE RELATIVISTIC a-u MODEL 

Tliere are several models to deal with the superdense matters, such as non-relativistic 
models, relativistic fieldtheoretical models [3] . To different models, the fracions of particles 
in the superdense matters are different, and then the bulk properties of superdense matters 
are different, that is, the EOS of them are different. Here the Relativistic a — u Model 
will be adopted [10]. The Lagrangian density of this model is 

L = Y.'^siilij.d^ + mB- QaBCr - gc,Blij.UJ^ - \gpBlnT.p^)ll)B + 
B ^ 

\{Q<yf - \mW - \f,.F^^ + \mlu:,u:^ - U{a) - 

\p,.-r + \mlp,.p^ + E Uil,.d^ - mi)^i (12) 

in which U{a) ^ aa + ^ca^ + ^da"^, F^^ = d^u^ - d^u^, p^^ = - is 
the field operator of Baryon B{B—n, p for the traditional neutron stars, B—n, p. A, 
E, S, A for the hyperon stars); ipi is the field operator of lepton l{l=e, p); and a, uj'^, 
^ are the field operators of a, ui, p meson respectively; Qo^b, duiB, QpB are the coupling 
constants between ex, a;, p meson and baryon B respectively In general, the coupling 
constants between a{u} or p) meson and neutron and proton are equal and decided by 
saturated property of nuclear matter (c/o-n = Qap) Qum — 9uip) and the symmetry energy of 
nuclear matter(g'o-p, 9'wp); ?7iB,mi,mi, (i = a,u!,p) are the mass of baryon, lepton, meson 
respectively; r is the isospin operator. To leptons, we assume they are free fermi gas. Prom 
this Laglangian density, we can obtain the EOS of the superdense matters, p — p{p), in 
which, p and p are the pressure and energy density of the superdense matters, respectively. 
In 1-loop approximation, the loop's contribution to the propagators of nucleon and a 
meson is considered, and the renormalization is used to renormalizing the divergent part 
of loop contribution. In the numerical calculation, we adopt the value of the parameter as 
following [10]: a = -2.1 x WMeV\^ c = 0.97 xM^, d = 1277, Qs = 6.73, = 8.59, M„ = 
938MeV, = 783Mel^, m„ = 550MeV, rUp — 770MeV, and the incompressibility of 



nuclear matter is 22AMeV, which is consistent with the experiment result [11,12]. From 
Fig.l, it is clear that the EOS of traditional neutron stars is stiffer than the EOS of 
hyperon stars. 

Fig.l 

IV. NUMERICAL RESULTS AND DISCUSSION 

In order to compare the numerical results with the observation, we present some typ- 
ical observation value here. As we know, only a few masses have been determined in 
observation from the more than thousand neutron stars. There are two typical observed 
mass data: 1.36 ± O.OSMq [13] for double neutron star binaries, and 1.87lo;i7-^Q(Vela 
X-1) [14] and 1.8 ± 0.4MQ(Cygnus X-2) [15] for X-ray binaries. 

The numerical results for the rotating traditional neutron stars and hyperon stars in 
the relativistic a — uj model are showed in the following figures and table. 

Figs. 2-3 show the changes of masses and radii between rotating neutron stars and 
non-rotating neutron stars as function of rotational period. At a given central density, it 
is easy to see that, as the neutron stars rotating slowly, the increment of the masses and 
radii will reduce sharply when the neutron stars' rotational periods are small than 1 ms. 
But as the periods become bigger than 1.6ms, which is the period of the fastest rotating 
pulsars in observation [3] , the change of the increment of the masses and radii caused by 
the rotation will become very week, and the increment will be no more than 2 percent. As 
the changes of the mass and radius to a real neutron star caused by the rotation are very 
small comparing with the total mass and radius, we can see that Hartle's approximate 
method is rational. 

Figs. 4-7 show the masses and radii of rotating and non-rotating traditional neutron 
stars and hyperon stars as function of central density at a given central angular velocity 
relative to the local inertial frame. From these figures, one can see that around the typical 
observational radii with value of 12km, the increments of the radii are bigger than that 
of other place. In table 1, some typical value in these figures are listed, combining fig.4, 
fig. 6 and the observational value, one can see that the double neutron star binaries are 
more like hyperon stars, while the neutron stars of X-ray binaries are more like traditional 
neutron stars. For the hyperon stars, its EOS is so soft that, even the period small than 
the smallest observational value, the total mass can't reach the observation value of X-ray 
binaries, about I.SMq. 
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Fig. 8 shows the period as a function of the central density at two giving central angu- 
lar velocities relative to the local inertial frame. From this figure, one can see that bigger 
central density and bigger central angular velocity relative to the local inertial frame cor- 
respond with a smaller period. In fig.8 one can find out that as the EOS of hyperon star 
is softer, at the same conditions, the period of hyperon star is appreciably bigger than 
the period of the traditional neutron star. Figs. 9-10 show the masses and radii of the 
traditional neutron stars and hyperon stars as a function of period. One can see that 
even the central density of the hyperon stars is bigger than that of the traditional neutron 
stars, the masses and radii of the hyperon stars are smaller than that of the traditional 
neutron stars. Another interesting result is that at a giving central density, when the 
period increases to a value which is bigger than the smallest observational period, 1.6ms, 
the change of the masses and the radii with the period is not obvious, which is consistent 
with the result of fig. 2. 

From these figures, one can see that at a giving period, the central density and the 
central angular velocity relative to the local inertial frame could be chosen freely, but 
if the masses and the period of a neutron star are given (by the observational value) at 
the same time, then the central density and the central angular velocity relative to the 
local inertial frame will be decided. As we know, there is another observational value: the 
radius of the neutron star, the decided central density and central angular are not always 
given the exactly radius. So one can say that if these three observational values: mass, 
radius and period are given to the same neutron star, there is only one special EOS could 
give a set of calculating values, which are fit like a glove to the set of observational value, 
that is, these three observational values could decide the EOS entirely. But the problem 
is that there is no neutron star with both mass and radius observationally determined up 
to the present. 



Fig. 2-Fig. 10 

Table 1. Rotating neutron star's property in the relativistic a — lu model* 
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TNS 


2.50 1.067 


11.98 


1.768 


12.18 


1.808 0.017 0.023 


1.488 




5.00 0.903 


10.21 


1.650 


12.94 


1.801 0.070 0.092 


0.791 


HS 


2.50 0.921 


11.80 


1.336 


12.04 


1.367 0.020 0.023 


1.686 




5.00 0.811 


10.87 


1.265 


12.87 


1.385 0.084 0.095 


0.873 



*In the table, TNS and HS denote the traditional neutron stars and the hyperon stars respectively; 
u)c is the angular velocity relative to the local inertial frame at the center, with a unit of (10^s~^); pc 
is the central density, with a unit of (lO^^fc^.m"^); i?o and Mq denote the radius and the mass of the 
non-rotating stars, R and M denote the radius and the mass of the rotating neutron stars, respectively; 
the unit of the mass is the solar mass(MQ), ^ and ^ are the fractional difference of the radius and 
mass between the rotating neutron stars and non-rotating neutron stars; P is the rotational period. 
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V. FIGURE CAPTIONS 

Fig. 1. Equation of state of traditional neutron stars(TNS) and hyperon stars(HS). 

Fig. 2. Increment of the mass and radius of TNS comparing to that of non- rotating 
TNS at the center density of pc — 3po, where po is the saturation density of nuclear 
matters. 

Fig. 3. Increment of the mass and radius of HS comparing to that of non-rotating HS 
at the center density of pc = 4po- 

Fig. 4. The total mass of the non-rotating and rotating TNS as a function of central 
density, where the mass is in a unit of solar masses. 

Fig. 5. The radius of the non-rotating and rotating TNS as a function of central 
density. 

Fig. 6. The total mass of the non- rotating and rotating HS as a function of central 
density, where the mass is in a unit of solar masses 

Fig. 7. The radius of the non-rotating and rotating HS as a function of central density. 

Fig. 8. Period as a function of the central density at two different angular velocities 
relative to the local inertial frame at the center. 

Fig. 9. The masses as a function of the period at two different central density for the 
TNS and HS. 

Fig. 10. The radii as a function of the period at two different central density for the 
TNS and HS. 
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